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Abstract 

We discuss the determination of the parameters of the pp and pp amplitudes for the 
description of scattering in the Coulomb interference region. We put enphasis on the 
possibility that the effective slope observed in the differential cross section is formed by 
different exponential slopes in the real and imaginary amplitudes (called Br and Bj). For 
this purpose we develop a more general treatment of the Coulomb phase. We analyse 
the differential cross section data in the range from 19 to 1800 GeV with four parameters 
(cr , p , Bj , Br ), and observe that we cannot obtain from the data a unique determination 
of the parameters. We investigate correlations in pairs of the four quantities, showing 
ranges leading to the smaller \ 2 values. 

In the specific case of pp scattering at 541 GeV, we investigate the measurements of 
event rate dN/dt at low \t\ [16] in terms of the Coulomb interference with exponentially 
decreasing nuclear amplitudes. The analysis allows a determination of the normalization 
factor connecting the event rate with the absolute cross section. 



I. INTRODUCTION 



The experiments in the high energy accelerators of Cern and Fermilab in the period 
from 1960 to 1990 collected data on the differential cross sections da/dt for the 
systems pp and pp at the center of mass energies a/s = 20 — 1800 GeV . Along almost 
50 years since the beginning of these studies of high energy hadronic scattering, many 
theoretical models were developed, but the experimental data stopped increasing in 
quantity or quality, as the accelerators were discontinued. The phenomenology and 
theoretical treatments of these systems are thus restricted in several aspects. The 
theoretical literature is enormous, now with increased interest due to the higher 
energy data that will come from LHC operation pp, and in the present work we do 
not analyse or compare the many dynamical and phenomenological efforts. 

To study of the dynamics that governs the processes, it is necessary to disentangle 
the squared moduli of complex quantities that represent the measured quantities 
in terms of the imaginary and real parts, which are intrinsically combined with the 
Coulomb contribution. 

In the region of small momentum transfers \t\, the intense Coulomb amplitude added 
to the nuclear interaction creates an interference that is observable in the \t\ distri- 
bution in da/dt. This Coulomb interference region of low \t\ values goes typically 
up to \t\ = 0.01 GeV 2 , but we show that the form of da/dt in general can actually 
be described in terms of simple exponential real and imaginary nuclear amplitudes 
well beyond this range. 

In previous analysis of the pp and pp data, the real and imaginary nuclear amplitudes 
were considered as having the same exponential dependence exp (Bt/2) , where B 
is the slope of the log plot of da/dt . This simplifying assumption is not adequate, 
according to dispersion relations j2] and according to the theorem of A. Martin [3] 
that says that the position of the zero of the real amplitude is close and approaches 
t = as the energy increases. Both results indicate that the slope of the real 
amplitude should be larger than that of the imaginary one, and in the present 
work we investigate the description of the Coulomb interference region allowing for 
different real and imaginary slopes. We review the scattering data in cases where 
this kind of information can be looked for. 

In Sec. II we review the expressions for the observable quantities in the forward 



2 



region, and obtain the expression of the intervening relative phase for the more 
general case of different slopes for the real and imaginary amplitudes. In Sec. Ill we 
analyse the differential cross sections for the energies in pp and pp scattering where 
data are more favorable. In Sec. IV we present some remarks and conclusions. 

II. LOW \t\ REGION AND COULOMB PHASE 
A. Description of scattering for small \t\ 

In elastic pp and pp collisions, the combined nuclear and coulomb amplitudes is 
written 

F c+N (s, t) = F c (s, t)e ia ^ + F N (s, t) , (2.1) 
where F° is the Coulomb part 

F° = (-/+) ~ F p 2 roton (2.2) 

with the proton electromagnetic form factor 

proton = (0.71/(0.71 + |t|)) 2 , (2.3) 

associated to a relative phase $, and F N is the strong interaction complexe ampli- 
tude 

F N (s,t) = F%(s,t)+iF?(s,t) . (2.4) 

The phase $ was initially studied by West and Yennie [5], and different evaluations 
have been worked out by several authors (6j [7J |8]. In the present work we extend 
these investigations considering the possibility of different slopes for the real and 
imaginary amplitudes. 

In the normalization that we use [I] the differential cross section is written 

^ = 7i\F c+N (s,t)\ 2 = 7r\F c (s,ty a ^ + Fg(s,t) +i F?(s,t)\ 2 . (2.5) 
For small angles we can approximate 

F N (s, t) « Fg{s, 0)e BRt/2 + tF?(s, 0)e B ' t/2 . (2.6) 

The slopes B R and Bj are usually treated as having equal values. In the present 
work we allow Br ^ Bj. 
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The parameter 



P 



F?(s,0) 



(2.7) 



the optical theorem 

a = 4:71 (0.389) Im Fj N (s, 0) , (2.8) 

and the slopes Br, Bj are used to parametrize the differential cross section for small 
\t\ . In these expressions, a is in milibarns and the amplitudes Fr, Fj are in GeV~ 2 



da 



-Bt 



For low \t\ , eq. (2.6) leads to the approximate form 

da 

~dt 

with 

B 



t=o 



dt 

P 2 B R + B r 



l + p 2 

as the usual slope observed in the data of da /dt. 



(2.9) 



(2.10) 



B. The Coulomb phase 

Here we derive an expression for the phase appropriate for cases with Br ^ Bj. 

The starting point is the expression for the phase obtained by West and Yennie [5] 

ro dt' 



*(*,*) = (-/+) 



In 



t\ 

Sj J-4p 2 \t' — t\ 



1 F N { Sl t') 



(2.H) 



F N (s,t) 

where the signs (— /+) are applied to the choices pp/pp respectively. The quantity 
p is the proton momentum im center of mass system, and at high energies 4p 2 w s . 

For small \t\, assuming that F N (s, if) keeps the same form for large \t'\ (this approx- 
imation should not have practical importance for the results), we have 
F N (s, f) F»(s, 0)e B ^'/ 2 + i F T N (s, 0)e^ i '/ 2 



F N (s, t) Fg{s, 0)e B ^/ 2 + i F T N (s, 0)e B ^/ 2 



c + i 



,B R {t'-t)/2 



+ 



,B z p-t)/2 



c + i 



where 



c = pe 



{Br-B^I/2 



(2.12) 



(2.13) 



The calculation is explained in detail in the appendix. The integrals that appear in 



the evaluation of eq. (2.11) are reduced to the form [8] 



1(B) 



r° dt' 


1 _ e B(t'-t)/2 


l-4p2 \t' - t\ 





(2.14) 



that is solved in terms of exponential integrals [U] as 



I(B) = Ei[*{4p l + t 



Ei[ - f ] + In [f (V + 1)} - In [ - ^] + 2 T . (2.15) 



The real part of the phase is then written 



<&(*,*) = (-/+) 



In 



J + c 2 + l 



c 2 J(S B ) + I(Bj) 



and this expression in introduced into eq. ( 2.5 ) . 



(2.16) 



III. ANALYSIS OF EXPERIMENTAL DATA 



With a in mb and t in GeV 2 the practical expession for da/dt in terms of the 
parameters a, p , Bj and Br is 



da 
~dt 



0.389 7T 



p a e 



B R t/2 



+ F C cos($) 



0.389 x 4vr 



a e 



B/t/2 2 

+ F c sin($) 



0.389 x 4vr 



where by $ we mean the real part given in eq. (]2.16|), written 

$ = (-/+) a In 

where 



1 1 z„ 



1 

1 + C 2 



c 2 I{Br) + I{Bi) 



with c given in Eq. (2.13) 



At high energies and small \t\ we simplify 



4p 2 + s 



and then the functional form of 1(B) is written 

r/ „. „ fBs\ „ ( Bt\ , [Bs 
HB) = El - )-EA-- Uln - 



hi 



fit 



+ 2 7 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



We have used fitting programs (Cern Minuit-PAW and Numerical Recipes) to obtain 
correlations for the four parameters ( a, p, Br, Bj ), for some values of energy where 
the data from CERN and Fermilab [TU] have more quality and quantity. Below we 
present these cases. 

It is important to remark that the results obtained in the fittings in general depend 
strongly on the set of data of low \t\ selected for the analysis of the Coulomb inter- 
ference region. This shows that the data accumulated in these experiments are not 



TABLE I: Forward scattering parameters found in the literature. At 1800 GeV the data 
values 1800( a ) and 1800^^ correspond to the experiments E710 and E741 in Fermilab 



(GeV) 


a (mb) 


P 


5(GeV~ 2 ) 


19.4 


38.98 ±0.04 


0.019 ±0.016 


11.74 ±0.04 


23.5 


38.94 ±0.17 


0.02 ±0.05 


11.80 ±0.30 


30.7 


40.14 ±0.17 


0.042 ±0.011 


12.20 ±0.30 


44.7 


41.79 ±0.16 


0.0620 ±0.011 


12.80 ±0.20 


52.8 


42.67 ±0.19 


0.078 ±0.010 


12.87 ±0.14 


62.5 


43.32 ±0.23 


0.095 ±0.011 


13.02 ±0.27 


541 


62.20 ± 1.5 


0.135 ±0.015 


15.52 ±0.07 


1800^ 


72.20 ±2.7 


0.140 ±0.069 


16.72 ± 0.44 


1800^ 


80.03 ±2.24 


0.15 


16.98 ±0.25 



detailed and regular enough to allow precise determination of the amplitudes in the 
forward direction. 

We stress that in this paper we do not intend to give new better values for pa- 
rameters. Instead, we show that the analysis of the data leads to rather ample 
possibilities. 

The values of forward scattering parameters for pp scattering at Fermilab and Cern 
ISR energies given in the standard literature are given in Table [!} The data at 
a/s = 19.4 GeV come from Fermilab , and that at 23.5 - 62.5 GeV are from the Cern 
ISR, with a review by Amaldi and Schubert [llj. 



Our determination of the parameters is described below. In Table [IT] we collect our 
results. It is remarkable that the obtained x 2 values are very small. 
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TABLE II: Forward scattering parameters obtained in our analysis. 



Vs (GeV) 


a(mb) 


P 


5/(GeV- 2 ) 


^(GeV- 2 ) 


P = B R /Bj 


x 2 


19.4 


40.379 ± 0.069 


0.019 (fixed) 


14.539 ± 0.262 


Bi , 25/ 


1 , 2 


1.299 


23.5 


39.821 ± 1.479 


0.0186 ±0.0137 


14.912 ± 9.246 


35.220 ± 177.57 


2.36 


0.2952 


30.7 


40.024 ±0.047 


0.027(fixed) 


11.784 ±0.239 


Bi , 2Bj 


1 , 2 


0.5361 


44.7 


41.839 ±0.291 


0.0543 ± 0.0037 


12.976 ± 0.631 


16.132 ± 15.509 


1.243 


0.6110 


52.8 


42.576 ±0.820 


0.0799 ± 0.0086 


13.414 ± 1.847 


14.113 ±33.684 


1.052 


0.1138 


62.5 


43.298 ±0.159 


0.0867 ±0.0034 


13.299 ± 0.358 


13.900 ± 10.007 


1.045 


0.5389 



A. pp scattering at yfs = 19.4 GeV 

Considering only Kuznetsov [12] and Schiz [TH] measurements, we have a total of 
69±134 points at 19.4 GeV. We have fittted the set of the first 61 points from 
Kuznetsov plus the 12 first ones from Schiz, covering the range 

0.00066 < |t| < 0.0395 . 

According to the information in the Durham Data Basis about this experiment, it 
is known that fittings that include the points of Kuznetsov lead no negative values 
of the parameter p. We then fix the value p = 0.019 (taken from Table [T]) and leave 
free the other parameters. The results are given in Table [IT] . The data and the 
solution of fitting are shown if Fig. [T] . It is remarkable that the choices Br = Bj 
and B R = 2 Bj lead to the same \ 2 - 

B. pp scattering at y/s — 23.542 GeV 

At 23.542 GeV there are 31 experimental points [H]. In order to obtain smaller 
values of x 2 , the fitting was made using the first 17 points, with \t\ in the interval 

0.00037 < |*| < 0.00395 , 

leading to the results given in Table [TT| The central value obtained for the ratio (3 is 

P = Br/Bj = 2.36 . 

The same value \ 2 = 0.2952 is obtained for any (3 = Br/Bj in the interval from 
1.22 to 4.91 . 
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FIG. 1: Data from Kuznetsov (61 points) |12j and Schiz (12 points) [13] are fitted to 
determine parameters at 19.4 GeV. We fix p = 0.019 (taken from Table |l| because fittings 
including Kuznetsov's measurements lead to negative value of p. 
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The plot of the 31 points together with the line obtained with fitting of 17 points 
are shown in Fig[2| Correlations between parameters p and j3 are shown in the RHS, 
with level curves of \ 2 ■ This plot is made fixing the values of p and (3, so that only 2 
parameters are free (this explains the relation 0.2952 = 0.25587 x (17 — 2)/ (17 — 4) 
in the values of x 2 ). The two different algorithms lead to distinct ranges in /?, 
indicating that the data have poor definition for this quantity. 

C. pp scattering at y/s = 30.632 GeV 

Although the data (32 points) of pp scattering at 30.632 GeV [TJ] look regular , we 
have difficulties to find values for the parameters, although x 2 comes out small. The 
value 0.042 given for p in Table|T]does not seem to be realistic, as our procedure leads 
to smaller values. With the value p = 0.027 fixed, the ratio (3 = B R /Bj can vary 
in a large interval, keeping the same \ 2 = 0.5361. The data, fitted with Coulomb 
interference expressions and exponential amplitudes, are shown in Fig. [3] and the 
parameters are given in Table [XTJ 
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FIG. 2: Data at 23.542 GeV . The first 17 points are fitted with expressions of Coulomb 
interference and exponential forms for the amplitudes. The graph in the RHS shows the 
regions of variation of p and (3 = Br/Bj that correspond to small values of x 2 while the 
quantities a and Bj are free. 




D. pp scattering at yfs = 44.699 GeV 

The 230 data points of the experiment at y/s = 44.699 GeV extending up to \t\ ~ 
7 GeV 2 are presented in the report by Amaldi and Schubert [llj. The forward part, 
with 40 points, up to \t\ ~ 0.02 GeV 2 , is very well fitted by the Coulomb interference 
formula, as shown in Fig. [4| with \ 2 = 0.6110. The parameters are given in Table 

M. 

Any value of the ratio (3 = Br/Bj in the range 

Br 

0.601 < -f- < 1.667 
Bi 

leads to the same value 0.6110 for x 2 - 111 Fig- p^we show the solution for low \t\ 
together with higher \t\ data, exhibiting the peculiar behaviour of the amplitudes 
deviating from the simple exponential dependence. In the plot with all 230 points , 
the dotted line shows the fitting obtained with the parametrization used in a previous 
work [12]. The parameters obtained in this case are a = 42.10 mb, p = 0.094, 
Bi = 12.22 GeV -2 , and B R = 24.43 GeV -2 . It is interesting that the description 

9 



FIG. 3: Data at 30.632 GeV . The 32 points are fitted with expressions of Coulomb 
interference and exponential forms for the amplitudes. 
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FIG. 4: Fitting of the 40 points with lowest \t\ at 44.699 GeV . 
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FIG. 5: The extended dashed line representing the fitting at very low \ t\ is plotted together 
with measurements at higher \t\, at y/s = 44.699 GeV. The dotted line shows a fitting of 
the whole set of points [15] . 




of the whole \t\ range leads to a definite indication for Br k 2 5j. Correlations 
between parameters are shown in Fig[6j with level curves of \ 2 ■ The two plots show 
respectively the correlations between p and (3 and between a and (3. of 

With fixed B R jB I = 2 we obtain X 2 = 0.6111 , with a = 41.819 ± 0.089 , p = 
0.0543 ± 0.0029 , Bj = 12.943 ± 0.406 . 

E. pp scattering at y/s = 52.806 GeV 

Fig. [7] shows the forward data at y/s = 52.806 GeV (34 points) (TJ] and our fit- 
ting with Coulomb interference expressions using the first 20 points. Correlations 
between parameters are shown in Fig. [8] , with level curves of \ 2 ■ The two plots 
are built fixing the two parameters in the axes while the other two parameters are 
found by fitting. The parameters obtained with 20 points in the interval 

0.00107 < \t\ < 0.02235 

are given in Table [Tf] . The ratio (3 = Br/Bj with any value in the range 

0.702 < Br/Bj < 1.214 
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FIG. 6: Ranges of values of parameters at y/s = 44.699 GeV that lead to small values of 
X 2 - The searches are made fixing the two values in the axes of the plots, while the other 
two parameters are free. 




FIG. 7: Data of pp forward scattering at 52.8 GeV fitted with Coulomb interference 
formulae. 
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FIG. 8: Correlations between pairs of parameters in pp scattering at 52.8 GeV. 




leads to the same value 0.1138 for x 2 ■ 

With fixed B R jBi = 2 we obtain x 2 = 0.1140 , with a = 42.543 ± 0.057 , p = 
0.0801 ± 0.0047 , Bi = 13.347 ± 0.214 . 



F. pp at ^ = 62.5 GeV 



Fig. [9] shows the data (138 points) [14] and result of our fitting with Coulomb 
interference expressions using the first 40 points of the set. Correlations between 



parameters are shown in Fig. 10 , with level curves of x ■ The parameter values 



obtained in fitting with 40 first points, in the \t\ interval 

0.00167 < \t\ < 0.051 



are given in Table |H} The same x 2 — 0.5389 is obtained with ratio (5 with any value 
in the interval 

0.78 < B R jB 1 < 1.25 . 

With fixed B R jB I = 2 we obtain X 2 = 0.5400 , with a = 43.264 ± 0.025 mb , 
p = 0.0869 ± 0.0017 , 5 7 = 13.230 ± 0.077GeV~ 2 . 
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FIG. 9: Data (138 points) of pp scattering at y/s = 62.5 GeV |14] shown together with 
our fitting of the first 40 points of the set with Coulomb interference expressions. 
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FIG. 10: Correlations between pairs of parameters ( p — (3 and p — a in pp scattering 
at 62.5 GeV. 
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FIG. 11: Event rate jTB] and differential cross section at 541 GeV 




G. Analysis of the event rate dN/dt at y/s = 541 GeV 

The lowest \t\ values reached in measurements of pp elastic scattering in the neigh- 
borhood of 540 GeV are reported with event rates [16] only. The cross-section values 
have not been determined otherwise in this low \t\ range, according to Durham HEP 
data basis. In the present work we use the Coulomb interference as a tool to find the 
normalization factor connecting event rate and differential cross-section. We find a 
very clear and precise connection, which is very important, as the event rate dN/dt 
has been measured with homogeneous accuracy, with many (99) points in a range 
of low \t\ values. Our simple procedure is to fit the dN/dt data with the expression 
for the Coulomb interference region, with and arbitrary multiplying normalization 
factor, and the four parameters that describe forward pp scattering. 



Fig. 11 shows the 99 points of the dN/dt measurements [16], of pp elastic scattering 
at 541 GeV and the values of the differential cross section after we determine the 
normalization by Coulomb interference. We have found the normalization factor 

— = — x 1 ( 3 5 ) 

dt dt 10.083 ±0.135 v ' ; 



In Fig. 12 we show the good agreement of da/dt at low \t\ , obtained from dN/dt 
by adjustment of the coulomb interference, with the data of G. Arnison et al [T7] in 
a \t\ range that partially superposes with the dN/dt event rate data. In the RHS of 
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FIG. 12: The normalized event rate data at 541 GeV re in good agreement with the da/dt 
data pH HE! [H] at the same energy and higher \t\ values. 




the same figure we include also the data of Bozzo et al. JT8J [19], including high \t\ 
values. 

As a test of consistency of this method of connection between event rate and absolute 
cross section we compare values of da/dt at the Cern/ISR energies, multiplied by 
an arbitraty normalization factor, with the Coulomb interference amplitudes. We 
find that this normalization factor is actually equal to one in all investigated cases. 

Observing locally the comparison between the prediction of da/dt obtained from 
the event rate dN/dt by adjustment to Coulomb interference equations and the 



data [HJ EE] we see that there is a discrepancy of a few percent. In Fig. 13 we show 



that the perfect matching is obtained with normalization factor 10.6, namely 

da fdN\ , 

— = [ — ) / 10.6 . 

dt V dt J ' 



In Fig. 14 we show the differential cross-section measured by Bernard et al. [20] 
and Abe et al. [21] which are compatible with each other, but do not match our 
normalized solution for dN/dt. To match, we have to introduce an arbitrary normal- 
ization factor 11.0 instead of 10.083 that we have determined. The ratio (3 = B R /Bj 
in the range 

0.898 < (3 = — ^ < 1.111 
B i 
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FIG. 13: Looking closely at the normalized cross section obtained with factor 10.083 
derived purely from Coulomb interference formulae, compared to Bozzo et al. data |18[ [T9] , 
we observe that there remains a displacement of a few percent. A more precise matching 
is obtained with a factor (dN/dt)/10.6 , as shown by the dashed line in the figure. 




leads to the same value 1.097 for x 2 - 

To show the influence of value of the ratio Br/Bj beyond these limits , we have 
fitted with fixed 

(3 = Br/Bj = 2 



obtaining parameter values shown in Table III 



We have also calculated with fixed normalization factor 10.6 , with results shown in 
the table. 



We see that the x 2 values do not vary strongly, showing that the data can be 
described, within errors, by scattering parameters in different ranges. 



Correlations among parameters are shown in Fig. [15] with drawing of level lines 
determined for low x 2 values. 

We have built a file with a continuous and non superposing set of points, being 
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FIG. 14: The measurements of Bernard et al. [20 and Abe et al. [2T] are compatible 
with each other, but do not match the other measurements and our conversion from the 
event rate, requiring an arbitrary matching factor 11.0 . 




TABLE III: Forward scattering parameters 



cr(mb) 


P 


^(GeV" 2 ) 


^(GeV^ 2 ) 


normalization 


x 2 


63.897 ±0.377 
63.651 ±0.909 


0.172 ±0.009 
0.1601 ±0.0168 


15.347 ±0.145 
15.163 ±0.118 


15.452 ±4.577 
2 Bi (fixed) 


10.083 ±0.135 
10.268 ± 0.358 


1.097 
1.115 


62.688 ± 0.037 
62.842 ±0.147 


0.1484 ±0.0025 
0.1459 ±0.0070 


15.375 ± 0.035 
15.237 ±0.085 


17.662 ± 1.441 
2 Bi (fixed) 


10.6 (fixed) 
10.6 (fixed) 


1.119 
1.126 



the first 59 points from normalized dN/dt (with normalization factor 10.6) and 121 
points from Bozzo et al. The 180 points form a regular \t\ distribution, which we fit 
with formulae from our previous work [15]. The results are shown in Fig. 16 The 
parameter values obtained in this fitting are a = 63.06±1.90 mb , p = 0.124±0.005, 
5/ = 13.88 ± 0.42 GeV" 2 , B R = 25.79 ± 0.77 GeV" 2 , with X 2 = 1-32. 
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FIG. 15: Correlations between the parameters p and (3 and between a and (3 that lead to 
low values of % 2 . In each case the other two parameters are let free while tables of x 2 are 
built with specified values for the two plotted parameters. Also the normalization factor 
10.083 is taken as fixed independently from PAW fitting program. 




H. pp scattering at -y/s = 1800 GeV 

Two independent experiments, both at Fermilab, measured da/dt for small \t\ at 
1800 GeV, with a discrepancy of about 10 percent in the normalization for the total 



cross section. This is shown in Fig. 17 Our determination of parameters (with fixed 
p = 0.14) is given in Table IV Although this is not a close determination, due to 
large variation bars, it is interesting that the lowest \ 2 are obtained with Br larger 
than Bj, for both experiments.. 



From the results shown in Table IV, with large differences in \ 2 values, we learn that 
the data from the E710 experiment are more compatible with the forward scattering 



basic expression (3.1) for da/dt than the CDF data. Another important observed 
result is that increasing simultaneously the values p and Br we may obtain the same, 
or lower, \ 2 '■ the real amplitude at \t\ = may be larger, but it decreases faster. It 
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FIG. 16: Data at 541 GeV selected and organized. The dotted line is a fit described in 
our previous work |15j . 




FIG. 17: Data and results at 1800 GeV . 
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is as if for a given dataset we should keep constant a product pe~^ BR / 2 ^ cfi , with a 
value \t\ e ff that is effective for the Coulomb interference region. For the E710 data 
we extract from the table that |t| e fj- ~ 0.075 ; for the CDF basis the value is rather 
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TABLE IV: Forward scattering parameters at 1800 GeV 



Experiment 


cr(mb) 


P 


B I (GeV~ 2 ) 


B R {GeV- 2 ) 


x 2 


E710 


72.748 ±0.186 


0.14 (fixed) 


16.297 ±0.039 


115.57 ± 164.20 


0.6020 


E710 


71.824 ±0.184 


0.14 (fixed) 


16.282 ±0.039 


Bi (fixed) 


0.6060 


E710 


72.651 ±0.186 


1.0 (fixed) 


16.284 ±0.039 


167.93 ±48.561 


0.5961 


CDF 


80.917 ±0.436 


0.14(fixed) 


16.988 ±0.087 


72.006 ±116.15 


1.771 


CDF 


79.982 ± 0.432 


0.14(fixed) 


16.981 ±0.087 


5/ (fixed) 


1.775 


CDF 


80.159 ±0.433 


1.0 (fixed) 


16.865 ±0.087 


85.730 ± 16.937 


1.705 



large \t\ eS ps 0.28 . We thus see that the values of p and B R at 1800 GeV are not 
determined from the existing data in this analysis. 

The Fermilab data on da/dt at 1.96 TeV are now launched in preliminary form, 
covering the t\ range from 0.26 to 1.30 GeV 2 . In Fig. 



we show these new data 



together with the old 1.8 TeV data. The LHS plot gives also the lines representing 
the fittings of the E710 data (dashed line) and of the CDF data (dotted line) as 



informed in Table |TVj with fixed p = 0.14. 

The RHS figure shows the fitting (solid line) of all these data (marked by the interpo- 
lating dashed line) using our previous parameterization expressions (15] . We obtain 
x 2 = 1.675 for the whole lot of data (51±26±27=104 points). Taking only the E710 
points together with the 1.96 TeV data (51±27=78 points) we obtain a beautiful 
representation with \ 2 = 0.45; taking the CDF points together with the new 1.96 
data (26+27=53 points) we obtain \ 2 — 0.74. These results, that depend on our 
particular parameterization, but are meaningful because they come from a search 
for continuity in the data, again seem to favor the t dependence and normalization 
of the E710 data. 



I. Amplitudes 



The lines representing da/dt in Figs. 16 and 18 are obtained from real and imaginary 



amplitudes that are shown in Fig. [T9J with their characteristic slopes and zeros . 

Some typical values of parameters obtained in this analysis of the full t dependence, 
that can be read off from the plots of the amplitudes, are given in Table fV) These 
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FIG. 18: Data on pp scattering at 1.8 and (preliminary) at 1.96 TeV. The lines in the LHS 
plot represent the fittings with the basic expression ( |3.1| ) : dashed for E710 and dotted 
line for CDF data. The RHS plot shows the fitting (solid line) of all the data (dashed 
line) with parameterization used in our previous work [15 . 




0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 

Itl (GeV 2 ) Itl (GeV 2 ) 



FIG. 19: Amplitudes of pp scattering at 541 and 1800 (including preliminary 1960) GeV 
normalized to 1 at \t\ = 0. The amplitudes are obtained with parameterization used in 
previous work |15| . Notice the linear scale in the LHS plot. The zeros of the real and 
imaginary amplitudes at 1800/1960 GeV occur respectively at about 0.05 and 0.68 GeV 2 
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TABLE V: Quantities extracted from the t dependence of the amplitudes at 541/546 and 
1800/1960 GeV. Remarks: (l)Abe+Amos+1.96 ; (2) Amos+1.96 ; (3) Abe+1.96 . 





a 


P 


Bi 


Br 


\t R \ 
l r o 1 


K\ 


x 2 


(GeV) 


(mb) 




(GeV~ 2 ) 


(GeV" 2 ) 


(GeV 2 ) 


(GeV 2 ) 




541/546 


63.05 


0.12 


13.88 


25.79 


0.16 


0.85 


1.32 


1800/1960 (1) 


73.98 


1.17 


15.50 


85.43 


0.05 


0.69 


1.68 


1800/1960 (2) 


73.95 


0.75 


15.47 


84.12 


0.05 


0.69 


0.45 


1800/1960 (3) 


88.49 


1.21 


17.94 


61.43 


0.05 


0.69 


0.74 



are possible, but not unique or necessarily correct, representations of the data. 



The real amplitude shows a large slope and falls rapidly to zero at a value \t$ \ that 
is expected to approach the origin in the form [15] 

following Martin's discussion about the first zero of the real amplitude [3]. Obviously, 
as |io| approaches the origin, Br grows with the energy. 

The numbers indicate much indetermination in this analysis of the old data. In 
particular, identification of the values of p and Br seems to be difficult. Dispersion 
relations must be used as a guide towards the disentanglement, in particular with 
application of the derivative dispersion relations for slopes 1 2] . We recall that the 
energy dependence in the measurements is crucial for the good use of dispersion 
relations, and we hope that LHC will produce diffractive data at several energies ? 



IV. REMARKS AND CONCLUSIONS 



We have performed a detailed analysis of the data on differential elastic cross-sections 
allowing for the freedom of different slopes for the real and imaginary amplitudes, 
namely Br ^ Bj, for the data points obtained in the ISR/SPS(Cern) and Tevatron 
(Fermilab) experiments during the years 1960-1990. The principal conclusion of the 
present analysis is that different values of the slopes, in particular the possibility 
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of B R > Bj in accordance with the expectations from Martin's theorem [2] and 
from dispersion relations [1], are perfectly consistent within the present errors of 
experimental data. 

Our investigation concerns the four quantities relevant for the elastic forward pro- 
cesses, namely, a,p,Bn and Bj. Studying the behavior of x 2 values near its min- 
imum and statistically equivalent parameter ranges, we observe that the available 
data from Cern and from Fermilab for small \t\ at the energies 20 — 2000 GeV are 
not sufficient for a precise determination of these four parameters. With real and 
imaginary amplitudes as independent quantities in the form a e' 3 *, the x 2 analysis 
clearly shows a strong correlation among the parameters, exhibiting a very large val- 
ley in x 2 surface in the parameter space. This in particular leaves large ambiguities 
in the determination for the weaker real part. 

It is interesting to note that detailed behavior in the x 2 surface in the parameter 
space depends very much on the chosen data sets. For example, the minimum of 
X 2 /degrees-of- freedom varies from 0.1 to 1.3 in different energies. On the other 
hand, the variation of x 2 f° r appreciable changes of some parameters is less than 0.1 
percent, in general indicating that their precise determination is not possible within 
the existing experimental situation. 

On the other hand, it is worthwhile to mention that, if we have sufficient number 
of data points in this region, we do not need the absolute normalization of the 
luminosity, since the Coulomb interference can determine correctly the absolute 
value of cross section as a free parameter, as was shown in this work for the 541 
GeV case. Of course, this method requires that the low \t\ data be very accurate. It 
is interesting to note that tests of this of Coulomb interference method, introducing 
a new free normalization parameter into the data set, are found to be compatible 
in both E-710 and E-741 Fermilab data at 1800 GeV, while the difference in the 
evaluation of the total cross section remains. 

As mentioned before, the quantitative analysis for the disentaglement of the real and 
imaginary parts of complex amplitude at small \t\ cannot be made with confidence 
with the data available up to now. Tests of quantities like the position of the zero of 
the real amplitude are not safe in these conditions. Thus model independent values 
for the four parameters cannot be obtained accurately. 

Usually, the models suited for pp dynamics aimed to cover an overall large \t\ region, 
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with correct description of dip and tail in da/dt, are not sensitive enough to the 
details of the behavior of scattering amplitudes for very small \t\ . In this region, the 
behavior of the amplitude maybe very sensitive to specific dynamical influences of 
the non-pertubative QCD dynamics. In particular, the forward scattering amplitude 
is directly related with the proton structure, and intimately related with the parton 
distribution function at small x and saturation problems. 

As a complimentary study, we also performed the analysis of scattering amplitude 
for a larger \t\ domain as proposed in [15]. We observe that such analysis, consistent 
with the general structure of the scattering amplitudes, such as positions of zeros and 
dips, clearly shows the necessity of the very distinct values of the slope parameters 
of real and imaginary amplitudes, corroborating with the results obtained here from 
the analysis in small \t\ regions. In terms of scattering amplitude, this different 
behaviors of real and imaginary part is crucial and will be fundamental to understand 
the mechanism of elastic scattering amplitude. 

Of course, the separation of complex amplitudes will never be a very easy task only 
from the measured scattering experiments, and it will be unavoidable to make use 
of other theoretical tools, such as dispersion relations. Particular attention must be 
given to the development and application of dispersion relations for slope parameters 
[1]. The knowledge of the energy dependence is crucial for the application of these 
tools in practice, and it is extremely interesting if an energy scan program is also 
included in the first phase of the LHC operation of p+p collisions. 

Although no new or better method has been introduced in the present work, our 
precise analysis has revealed the existence of correlations and related uncertainties 
of the behavior of the scattering amplitude at small \t\. For quantitative determi- 
nations, more precise and numerous data points are necessary. 

We expect to have a different situation in the future experiments from RHIC and 
LHC, with much better statistics and accuracy in the measurements of scattering 
data in the Coulomb interference region, together with a systematic energy scan 
program. 
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V. APPENDIX: THE CALCULATION OF THE COULOMB PHASE 



Here we give details of the evaluation of the phase of West and Yennie given by eq. 



(2.11 ) in the case where we let Br ^ Bj . 



After eq. (2.12) we define 



G 



c c(c — 1} 



R 



Gi 



c + i c 2 + 1 
1 (c-i) 



c + i c 2 + 1 



o dt' 



and the integral in eq.(2.11) is written 

F N (s,t'Y 
F N (s,t). 

1 _ e B R Xt'-t)/2 



1 - 



and then 



4p 2 \t' - t 
J-4p 2 \t' — t\ 



*(*,t) = ("/ + ) 



+ iGj 



dt' 

-4p 2 \f - t\ 



1 _ e Bi(t'-t)/2 



H--) + g r [ 

s J- 



+ iGj J 



o dt' 



-4p 2 \t' - t\ 

We note that both integrals are of the form 

r° dt' 
1(B) ' 



o dt' 

4p 2 \V - t 
1 _ e Bi(t'-t)/2 



1 _ e B R (t'-t)/2 



1 - e 



B(t'-t)/2 



(5.1) 



, (5.2) 



(5.3) 



(5.4) 



-4p2 \f - t\ 

which has been studied by V. Kundrat and M. Lokajicek [8]. With x = f — t and 
y = Bx/2, we have 



KB) 



o dx 
-4p 2 I a; I 



o dy 

-B(4p2 +t )/2 \y\ 
B(4p2 +4 )/ 2 

" \y\ l 



1 - e Bx ' 2 

l-e y 



+ 



1 - e 



-Bt/2 dy 

B(4p2 +t )/2 \y\ 
Bt/2 dy 

fi \y\ l 

Bt/2 d y 

\y\ 



l-e v 
l-e y 
e y - 1 



(5.5) 
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These expressions can be written in terms of exponential integrals, as can be seen 
in the Handbook of Mathematical Functions of M. Abramowitz and L.A. Stegun [§] 

as 



B(4p 2 +i)/2 d y 

\y\ 



1 - e 



V + t 



+ In 



4p 2 +t 



and 



e y - 1 



-Bt/2 rfy 

'o 1 2/ 1 

where 7 = 0.5772 is the Euler constant 



K(-ft)-ln(-ft)- 7 



+ 7 (5.6) 



(5.7) 



Thus the integrals that appear in eq. (5.2) have the functional form 
B 



Bt, 



1(B) = E 1 [- (4p 2 + f J] - Ei[ - ^] + In [- (V + t)] - In [ 
and the phase can be written 

t) = (-/+) In ( - -^j + G R I{B r ) + tGj I(Bj) 

or 



Bt 
~2 



] + 2 7 , (5.8) 



9{a,t) = (-/+) 





1 








+ c 2 + 1 



c 2 + l 



(5.9) 



1(5,) - J(S fl ) 

(5.10) 



The real part of the phase is taken into eq. ( 2.5 
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